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Abstract —Utilizing highly synchronized measnrements from 
synchrophasors, dynamic state estimation (DSE) can he applied 
for real-time monitoring of smart grids. Concurrent DSE stud¬ 
ies for power systems are intolerant to unknown inputs and 
potential attack vectors — a research gap we will fill in this 
work. Particularly, we (a) present an overview of concurrent 
estimation techniques, highlighting key deficiencies, (b) develop 
DSE methods based on cubature Kalman filter and dynamic 
observers, (c) rigorously examine the strengths and weaknesses 
of the proposed methods under attack-vectors and unknown 
inputs, and (d) provide comprehensive recommendations for 
DSE. Numerical results and in-depth remarks are also presented. 

Index Terms —Cubature Kalman filter, cyber attack, dynamic 
state estimation, observer, PMU, unknown input. 

I. Introduction and Motivation 

S TATE estimation is an important application of the energy 
management system. However, most widely studied static 
state estimation methods Q, 12 , a, a, Q cannot effectively 
capture the dynamics of power systems due to its dependency 
on slow update rates of supervisory control and data acqui¬ 
sition (SCADA) systems. In contrast, dynamic state estima¬ 
tion (DSE) enabled by phasor measurement units (PMUs), 
devices with high update rates and high global positioning 
system (GPS) synchronization accuracy, can provide accurate 
dynamic states of the system, and will thus play a critical role 
in achieving real-time wide-area monitoring, protection, and 
control in smart grids. The main objective of North American 
SynchroPhasor Initiative (NASPI) is to improve the visibility 
and reliability of power systems by wide area measurement 
and controlH Subsequently, efficient DSE methods can be 
utilized to improve the visibility of power networks. 

The main goal of this paper is to present DSE alternatives 
that address major limitations of current DSE methods such as 
intolerance to unknown inputs, attack-vectors, and parametric 
uncertainty. To achieve that, we study DSE methodologies 
for power systems, utilizing nonlinear, dynamic observers and 
recently developed, efficient Kalman filters. Many of the recent 
research efforts on DSE heavily focus on Kalman filter-based 

A. F. Taha is a research intern with the Energy Systems Division, Argonne 
National Laboratory, Lemont IL 60439 USA in 2015 and a Ph.D. candidate at 
the School of Electrical and Computer Engineering, Purdue University, West 
Lafayette, IN 47907 USA (e-mail: tahaa@purdue.edu). 

J. Qi and J. Wang are with the Energy Systems Division, Argonne 
National Laboratory, Lemont, IL 60439 USA (e-mail: jqi@anl.gov; jian¬ 
hui . wang @ anl. gov). 

J. H. Panchal is with the School of Mechanical Engineering, Purdue 
University, West Lafayette, IN 47907 USA (e-mail: panchal@purdue.edu). 
'available at https://www.naspi.org/. 


methods for reasonable state estimation. We show in this paper 
that this approach has its own limitations, in comparison with 
observers for nonlinear dynamical systems. 

Consequently, our objective is not to develop theory for state 
estimators or observers for generic dynamic systems, but to (a) 
survey the literature for current DSE methods, (b) discuss the 
limitations of the current status quo of DSE, (c) present and 
apply better alternatives, and (d) discuss the shortcomings and 
strengths of the proposed solution. 

II. Literature Review and Paper Contributions 

Here, we present an overview of DSE methods, introduce 
dynamic observers, and highlight our main contributions. 

A. Kalman Filters for Dynamic State Estimation 

Unlike many estimation methods that are either computa¬ 
tionally unmanageable or require special assumptions about 
the form of the process and observation models, Kalman filter 
(KE) only utilizes the first two moments of the state (mean and 
covariance) in its update rule 0, 0. It consists of two steps: 
in prediction step, the filter propagates the estimate from last 
time step to current time step; in update step, the filter updates 
the estimate using collected measurements. 

KE was initially developed for linear systems. However, 
for DSE the system equations and outputs have strong non¬ 
linearity. Thus variants of KE which can deal with nonlin¬ 
ear systems have been introduced. Eor example, DSE has 
been implemented by extended Kalman filter (EKE) IHl, ID, 
unscented Kalman filter (UKE) m, mi, m, square-root 
unscented Kalman filter m, Qa, extended particle filter ua, 
and ensemble Kalman filter ina. The reader is referred to uni 
for a survey of DSE techniques. 

B. Observers for Dynamical Systems and Power Networks 

Dynamic observers have been thoroughly investigated for 
different classes of systems. To mention a few, they have been 
developed for: linear time-invariant (LTI) systems, nonlinear 
time-invariant (NLTI) systems, LTI and NLTI systems with un¬ 
known inputs, sensor and actuator faults, stochastic dynamical 
systems, and linear and nonlinear hybrid systems. 

Most observers utilize the plant’s outputs and inputs to 
generate real-time estimates of the plant states, unknown 
inputs, and sensor faults. The cornerstone is the innovation 
function — sometimes a simple gain matrix designed to nullify 
the effect of unknown inputs and faults. Linear and nonlinear 
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functional observers, sliding-mode observers, unknown input 
observers, and observers for fault detection and isolation are 
all examples on developed observers for different classes of 
systems, under different assumptions. 

In comparison with KF techniques, observers have not 
been fully utilized in DSE. However, they inherently possess 
the computational, technical, and theoretical capabilities to 
perform good estimation of the grid’s physical status. 

C. Differences between Kalman Filters and Observers 

Observers are typically used for deterministic systems, 
whereas KFs are used for stochastic dynamical systems. If 
statistical information on process and measurement noise is 
available, KFs can be utilized by assuming a certain distribu¬ 
tion of noise. However, this assumption can be very strict as 
quantifying distributions of noise is challenging. 

While some deterministic observers may be robust and toler¬ 
ant against process and measurement noise, KFs are inherently 
designed to deal with model approximations and measurement 
noise. As for implementation, observers are simpler than KFs. 
For observers, matrix gains are computed offline to guarantee 
the asymptotic stability of the estimation error. 


D. Paper Contributions 


The contributions of this paper are summarized as follows. 

1 ) Surveying the most widely used DSE methods and dis¬ 
cussing their shortcomings; 

2 ) Extending existing DSE methods by introducing Cubature 
Kalman hlter (CKE), which directly uses nonlinear model 
as UKE but has better numerical stability; 

3 ) Presenting two dynamic observers for nonlinear systems 
and applying them to power systems; and 

4 ) Exhaustively studying the experimental and theoretical 
limitations and strengths of the proposed estimation meth¬ 
ods under different adverse scenarios. 


The paper is organized as follows. In Section III we repro¬ 
duce the nonlinear dynamics of the single-machine inhnite bus 
system. The CKE and two dynamic observers are introduced 
in Sections |rV| and |V] Then, numerical results are presented 
in Section |VI| Einally, insightful remarks, recommendations. 


and conclusions are presented in Sections VII and VIII 


HI. Single-Machine Ineinite-Bus System Model 


Here, we introduce the dynamics of the single-machine 
inhnite bus (SMIB) power system, which has been used as 
a benchmark in DSE literature, such as 13, ifTSl , lfT3 . Let 
X = [a;i X2 xz X4j = [6 uj e'^ e^] be the internal state 
of the system, where S is the rotor angle (rad); oj is the 
rotor speed (rad/s); and are the transient voltage along 

q and d axes. Also, let u = [mi 11,2]^ = [Tm be 

the input to the SMIB system, where is the mechanical 
torque (p.u.) and Ef^ is the internal held voltage (p.u.). 
The outputs are the real and imaginary part of the voltage 
phasor and current phasor obtained from PMUs installed at 
the terminal bus of the generator. We dehne the output as: 


y = [yi y2 ys 2/4] = [sr ej Ir ij] . The fast sub¬ 

transient dynamics and saturation effects are ignored and the 
generator model is described by the fourth-order differential 
equations in the local d-q reference frame, as follows: 


a::i = a:2 - ujq 

±2 = ^ {ui - Te - ^(X2 - Wo)) 
^3 = 7 ^ ( m2 - 3^3 - {Xd - x'^)id) 

X4= i^{-X4 + {Xq- X'q)iq) 


( 1 ) 


where is the electric air-gap torque; wq is the rated value 
of rotor speed; H and Kr are the inertia constant and the 
damping factor; iq and id are stator currents; T'g and T^g 
are open-circuit time constants; Xq and Xd are synchronous 
reactance; x'q and x'^ are transient reactance at the q and d 
axes, respectively. The iq, id, and Te variables in ([^, and the 
outputs 2/3 and 2/4 can be written as functions of x and u: 
'^R = Xi sin a;i + x^ cos a;i 
'k/ = xa sin a;i — X4 cos xi 

It = y(4- - e'qj) 

y3 = 'iR = Re(/t) 

y4 = ij = Im(/t) 

■ , 

Iq = [II Sin Xi + IR cos Xi ) 

ON 

Sb . . X 

Id= Or sin Xi - iicosxi) 
bN 

€q — X3 X^ id 
ed = X4 + x'g iq 




Sb 


^aio H" 


Gdid) 


Q 

ON 

where ^ is the generator voltage source, is the transient 
voltage along q axis of the inhnite bus, y is the line admittance 
of the reduced network which only has one line connecting the 
generator and the inhnite bus, and Sb and Sn are the system 
base MVA and generator base MVA, respectively. The outputs 
cr and ej can also be obtained as function of x and u: 

2/1 = Cr = X4 sin xi + Cq cos Xi (2) 

2/2 = 6/ = X3 sin xi — Cd cos xi. ( 3 ) 

Without unknown inputs, actuator faults, or disturbances, the 
state dynamics can be written in compact form as: 
x{t) = f{x, u) = Ax(t) + Buff) + 4 >{x, u), 

where 


A = 


0 
0 


1 

2H 
0 0 


0 0 


0 
0 

1 

' T' 

0 - 


4>{x,u) = 


qO 



' 0 

0 ■ 


Wo 

0 

1 
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2H 
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T ' 

^dO 



0 
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Wo 
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—Wo 


{-Tffx,u) + Ko) 


Xd - 

T ' 

^ do, 


T ' 
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id{x, u) 
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We utilize the above form of the power system dynamics in 
the estimation problem in the next sections. We also consider 
unknown inputs in the system dynamics and attack-vectors 
against output measurements. 

Remark 1. Albeit only the SMIB system is presented, the 
estimation methods in this paper can be applied to other 
systems. More details on the model for multi-machine system 
can be found in /EH?. We use this simple system as it is the 
benchmark in many recent DSE methods as in m and nss- 

IV. Kalman Filter Techniques 

As mentioned above, the current status quo of DSE is 
based on different variants of KF that can deal with nonlinear 
models, among which EKE, UKF, and especially CKF are 
discussed in this section. 

A. Extended Kalman Filter 

Although EKE maintains the elegant and computationally 
efficient recursive update form of KF, it works well only 
in a ‘mild’ nonlinear environment, owing to the first-order 
Taylor series approximation for nonlinear functions ll 20 l . It 
is sub-optimal and can easily lead to divergence. Also, the 
linearization can be applied only if the Jacobian matrix exists 
and calculating Jacobian matrices can be difficult and error- 
prone. For DSE, EKE has been discussed in IHl and Q- 

B. Unscented Kalman filter 

The unscented transformation (UT) ED is developed to 
address the deficiencies of linearization by providing a more 
direct and explicit mechanism for transforming mean and 
covariance information. Based on UT, Julier et al. El, Q 
propose the UKF as a derivative-free alternative to EKE The 
Gaussian distribution is represented by a set of deterministi¬ 
cally chosen sample points called sigma points. The UKF has 
been applied to power system DSE, for which no linearization 
or calculation of Jacobian matrices is needed Cni, ED- 

However, for the sigma-points, the stem at the center (the 
mean) is highly significant as it carries more weight which 
is usually negative for high-dimensional systems. Therefore, 
in terms of the numerical instability, the UKF is supposed to 
encounter troubles when used in high-dimensional problems. 

C. Cubature Kalman filter 


A nonlinear system (without unknown inputs or attack- 
vectors) can be written in discrete-time form as 

— f — — T Qk—i 

r/fc = h{xk,Uk) + Vk 

where Xk S K", Uk S K", and G are states, inputs, 
and observed measurements at time step fc; the estimated mean 
and estimated covariance of the estimation error are m and 
P; / and h are vectors consisting of nonlinear state transition 
functions and measurement functions; q^._i ~ N{ 0 ,Qi._i) is 
the Gaussian process noise at time step A: — 1 ; ^ -^( 0 , Rk) 

is the Gaussian measurement noise at time step k; and Qk_i 
and Rk are covariance matrices of Qk-i and r^. 

The procedure of CKF consists of a prediction step and an 
update step, which are summarized in Algorithms and 




Algorithm 1 CKF Algorithm: Prediction Steps 
I: draw cubature points from the intersections of the n 
dimensional unit sphere and the Cartesian axes. 

A/nCj, i = 1, • • • , n 
- Ci-n, i = n 1 , - ■ ■ , 2 n 
where Ci is a vector with dimension n, whose ith element 
is one and the other elements are zero, 
propagate the cubature points. The matrix square root is 
the lower triangular cholesky factor 

Xi,k-l\k-l = Pk-l\k-l $,i + 'mk-l\k-l- 

evaluate the cubature points with dynamic model function 


-1 — f{^i,k-i\k-i)- 


i,k\k 

4 : estimate the predicted state mean 

^ 2n 

■m^k-i = 7^^^*,k\k-i- 


2=1 


5 : estimate the predicted error covariance 


- 2n 

i —1 

- + Qk-i- 


V. Utilizing Dynamic Observers for DSE 

Here, we compactly present two observers that can be 
used for DSE. First, we develop the system dynamics under 
state disturbances, actuator faults, unknown inputs and attack 
vectors. Then, the two observers are succinctly presented. 


Both EKE and UKF can suffer from the curse of dimension¬ 
ality and the effect of dimensionality may become detrimental 
in high-dimensional state-space models especially when there 
are high degree of nonlinearities in the equations that describe 
the state-space model EOl . EH- Making use of the spherical- 
radial cubature rule, Arasaratnam et al. EOl propose CKF, 
which possesses an important virtue of mathematical rigor 
rooted in the third-degree spherical-radial cubature rule for 
numerically computing Gaussian-weighted integrals, better 
addressing the numerical instability problem in UKF. 


A. Dynamics Under Unknown Inputs and Attack Vectors 

Under disturbances, actuator faults and attack vectors, the 
dynamics presented in the previous section can be written as 
= fix, U) + BdUdit) -f Bafait) ( 4 ) 

y(t) = Cx{t) ( 5 ) 
where Bd is the disturbance matrix; Ba represents the actuator 
fault’s matrix; C denotes the output matrix (see Remark 1^1 
and v[f) exemplifies a potential attack vector against PMU 
measurements (see Remarks [H and 1^1. 
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Algorithm 2 CKF Algorithm; Update Steps 
1 : draw cubature points from the intersections of the n 
dimensional unit sphere and the Cartesian axes. 

2 : propagate the cubature points 

^i,k\k-l = Pk\k-l + f^k\k-l- 

3: evaluate cubature points using measurement function 


i,k\k — l ■ 

4: estimate the predicted measurement 

^ 2n 

yk\k — l ^ ^ ^ 1 ■ 


2n 

i=l 

5: estimate the innovation covariance matrix 

2n 


^k\k—l ^ ^, ^iMk — l'^i,k\k—l 


2n 

~ ifklk-iyklk-l + Rk- 
6 : estimate the cross-covariance matrix 

^ 2n 

Pxy,k\k-1 ^i,k-l\k-l^i,k\k-l 


2=1 


-T 

'^k\k — iyk\k—l ■ 

7 : calculate the Kalman gain 

I^k = Pxy,k\k-1^ k,k-l- 

8 : estimate the updated state 


TWfcife = mk\k-i + KkiVk - iJk)- 
9: estimate the updated error covariance 

Pk\k ~ Pk\k—1 ^kPyy,k\k — l^k • 


Remark 2. For simplicity, the output function is assumed to 
be linear in terms of the internal state of the system. Numerical 
results indicate that the differences between the outputs from 
nonlinear and linear output equations, given the same initial 
conditions and control input, of the SMIB system are minor. 

Remark 3. The National Electric Sector Cybersecurity Orga¬ 
nization Resource (NESCOR) developed cyber-security fail¬ 
ure scenarios with corresponding impact analyses l\24V . The 
report classifies potential failure scenarios into different 
classes, including wide area monitoring, protection, and con¬ 
trol (WAMPAC) — the paper’s main focus. The following 
WAMPAC failure scenario motivate the research presented in 
this paper: Measurement Data (from PMUs) Compromised 
due to PDC Authentication Compromise l\24]l . Hence, the 
addition of attack-vectors against PMU measurements depicts 
the aforementioned WAMPAC failure scenario. 

Remark 4. The attack vector, v(t), is different from typical 
measurement noise. While measurement noise vectors are often 
assumed to follow a certain distribution with small magni¬ 
tudes, the assumed attack vector here follows no statistical 
distribution, as demonstrated in the numerical results section. 
Measurement and process noise are also added. 

The unknown inputs to the given system are Ud (represent¬ 
ing the unknown plant disturbances) and (actuator faults). 
We assume that A, B, Bd, B^, and C are all known. For 


simplicity, we can combine the unknown inputs Ud and f ^ 
into one unknown input quantity, d{f), defined as follows 
d{f) = \ul(t) G The matrices Bd and 

are combined into D G and the plant dynamics can be 

rewritten as: 

x{t) = Ax(t) -f Bu(f) -\- 4>{x,u) -\- Dd{f) (6) 

y{t) = Cx{t) Fv(t). (7) 

Remark 5. The exciter output voltage Ef^ may not available 
for state-estimation. In that case Ef^ can be augmented to 
the unknown input vector and hence, matrices Bi and D 
will have different dimensions and values. Ghahremani and 
Kamwa apply EKE with unknown inputs to generate state 
estimates where the unknown input is Ef^ 4251? . 

Assumption 1. The nonlinear component in SMIB dynamics, 
4>{x, u), satisfies locally Lipschitz condition in a region D that 
includes the origin with respect to x, uniformly in u. Precisely, 
the nonlinear function satisfies the following condition: 

\\cl){x,u) - (t){z,u)\\ < 7 ||a; - z||, 
for all X and z in D, where y > 0. If H is equivalent to K”, 
the function (f>(x,u) is said to be globally Lipschitz, and it 
then satisfies the one-sided Lipschitz condition: 

{(f){x,u) - (t>{z,u),x- z) < p||a; - 2 ;f, 
where p G M. 

Remark 6. Although we define v{t) to be an attack vector, this 
definition is not restricting. The unknown quantity v(t) reflects 
possible measurement noise or falsely reported measurements. 
For example. It has been reported that PMUs from multiple 
vendors might produce conflicting measurements, as high¬ 
lighted in a North American Synchrophasor Initiative (NASPI) 
report (ESI?. Furthermore, the term Dd{t) models a general 
class of unknown inputs such as: nonlinearities, modeling un¬ 
certainties, noise, parameter variations, unmeasurable system 
inputs, model reduction errors and actuator faults sm, (EH?- 


B. Observer Design Using Linear Matrix Inequalities (EMI) 


Here, we introduce an observer design for a specific class of 
nonlinear systems with unknown inputs, which is based on the 
methods in 1^ by Chen and Saif. The observer design only 
assumes that D is full-column rank and that the nonlinear 
function in the state dynamics satisfies Assumption The 
observer design problem is then formulated as an LMI. 

Define & to be a real-time estimate of the real state of the 
system and 2 ; to be an intermediate state of the estimator. The 
dynamics of the observer developed in 1291 can be written as 
z{t) = Nz{f) -\- Gu{t) -\- Ly{t) -\- Mcj)(x, u) ( 8 ) 
x{t) = z(t) - Ey(t) (9) 

where matrices E and iC are design parameters and A?, G, L 
and M are obtained from the matrix equalities while guar¬ 
anteeing the convergence of the state estimates. Algorithm 1 
is applied to compute the aforementioned matrices. Given 
A, B, C, D and 7 , the above algorithm finds a solution to the 
observer dynamics by solving the formulated LMI ([T0|). The 
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Algorithm 3 Observer with Unknown Input Design Algorithm 
1 : compute matrices U,V,A and B: 

U^-D{CD)^ V = I - {CD){CD)^ 
A={I + UC)A B = VC A 
2 : find matrices Y, K and a symmetric positive definite 
matrix P that are a solution for this LMI: 

T’S'ii ^12 




12 ^ 2 n 


< 0 


( 10 ) 


where 




^11 = A P + PA 


B^Y^YB-C^ -kC+-fI, 

^12 = Vi {Pii + kC) + Yiyc)) 

obtain matrices Y and K and the observer parameters 
AT, G, L and M: 

Y = P-^Y , K = p-^k 
E = U + YV, M = I + EC 
N = MA- KC, G = MB, (11) 

L = K{I + CE) - MAE 


designed matrix parameters o guarantees that the observer 
dynamics yields asymptotically stable estimation error. 


Utilizing the solution to the LMI ( |T3] l, the state estimates 
generated from ( |T^ are guaranteed to converge to the actual 
values of the states. 

VI. Numerical Results 

In this section, we test the proposed estimation methodolo¬ 
gies on the SMIB system under different conditions. 

A. System Description and Parameters 

As shown in Fig. the SMIB system is extracted from PST 
toolbox 1321. We apply a three-phase fault at bus 3 of line 3—2 
(one of the parallel lines) to generate dynamic response. The 
fault is cleared at near and remote ends after 0.05s and 0.1s. 
DSE is performed for the post-contingency system, for which 
the parameters can be found in for the SMIB system. 


13 2 



0.1 

0.3999 



0.3999 

Infinite 



Bus 



h 


Generator 


Fig. 1. Single machine infinite bus power system. 


C. Observer Design for Nonlinear Systems 


Here, we present a recently developed observer in ll^ that 
can be applied for DSE. The examined observer assumes the 
nonlinear function 4){x,u) satisfies Assumption 1. In addition 
to this assumption, we assume that the nonlinear function is 
quadratically inner-bounded. Precisely, 

u) - cj){z, u)V {ct){x, u) - cj){z, u)) <V\x-zf 


{ct){x, u) - 4){z, u), X - z), 

where p and p are real numbers. Similar results related to the 
dynamics of multi-machine power systems established a sim¬ 
ilar quadratic bound on the nonlinear component (see. ED). 
The dynamics of the observer can be written as follows: 
x{t) = Ax[t) + Buit) + 4>{x, u) -f L{y{t) — Cx{t)), (12) 
where L is a matrix gain determined as follows. Eirst, given 
the Lipschitz constant^ p, p, and p, the following linear 
matrix inequality is solved for positive constants ei, e 2 and 
cr and a symmetric positive semi-definite matrix P, 

PA + {eip + e2p)In 
-aC^C 


A^P 




pe2 - €i 
_ 


< 0 . 


(13) 


Then, the observer gain-matrix L is given by: 

L^^P-^C^. 


^Those constants can be determined by evaluating the bounds on the 
nonlinear function (f>(x,u) and adding a possible bound on the unknown 
inputs and disturbances. Determining those constants affects the design of 
the observer, and hence it is advised to choose conservative bound-constants 
on the nonlinear function. The nonlinearities present in the SMIB system 
are bounded (e.g., sines and cosines of angles, multiplications of bounded 
quantities such as voltages and currents) 


This system has two known inputs (Tm and Ef^) and four 
outputs and ij). The two known inputs are kept 

unchanged. Without loss of generality, a system input (such 
as Ef^) can be augmented to the unknown input vector and 
considered as an unknown quantity (see Remark]^, in case it 
is not available. The SMIB’s and observer’s initial conditions 
are set as; 

a;(0) = [0.4233 377.3780 1.0277 0.1190]^, 

x(0) = [0 Wo 0 0]^ , 

where a:(0) is the post-contingency system state. The state- 
space matrices A, B, and C in 0 and 0 are omitted 
for brevity. The unknown input coefficients D in 0 are 
randomly chosen and the unknown input vector d{t) = 
[0.01cos(<) O.Olcos(f) O.Olcos(f)]^ (we consider three un¬ 
known inputs with the same magnitude). 

To test the robustness of CKE and the observer, a time- 
varying attack-vector v{t) against the four measurements is 
added to the measured output: 

Vi{t) = 0.2 sin(f) 


V 2 {t) = 


0.5sin(2f) -I- 

0 , 


5 < f < 6 
otherwise 


(14) 


vsit) = 0.3cos(f) 


V4{t) = 0 . 2 . 


This attack vector is manually chosen, showing different 
scenarios when the measured signals can be altered without a 
clear, known disturbance or attack signal (see Remark]^. 

Gaussian process and measurement noise are added. The 
process noise covariance is a diagonal matrix, whose diagonal 
entries are the square of 10% of the largest state changes, as in 
Ha, and therefore Q = diag(10-6,10-^ 1.47x10-1°, 1.17x 
10“®). The measurement noise covariance is i? = 10 -'iJ4 
where I 4 is an identity matrix with dimension 4. 
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Fig. 2. Comparison between the actual states and estimated ones for the 
SMIB system utilizing CKF without unknown inputs. Green and red dotted 
lines denote the case without or with cyber attack, respectively. 


Fig. 4. Attack-vector estimation for the SMIB system utilizing CKF. Green 
and red dotted lines denote the case with or without unknown inputs, 
respectively. Attack-vector is added for the two cases. 
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Fig. 3. Comparison between the actual states and estimated ones for the 
SMIB system utilizing CKF with unknown inputs. Green and red dotted lines 
denote the case without or with cyber attack, respectively. 


B. CKF Estimation Results 

The initial estimation error covariance matrix is chosen as 
a diagonal matrix whose diagonal entries are the square of 
1 % of the initial estimated states. If an initial estimated state 
is zero, the corresponding entry is set to be 0.01^. Therefore, 
Pq = 14.2122,10“"^, 10“'^). The estimation of the 

states is shown in Figs. and It is seen that in the absence 
of attack-vectors and unknown inputs, the estimates from CKF 
can converge to real states. However, when there are unknown 
inputs and/or cyber attack, CKF-based estimator fails to track 
the real states, indicating that CKF is not good at dealing with 
unknown inputs or cyber attack although it can filter noise. 

The estimation results of attack vectors for the cases with 
or without unknown inputs are shown in Fig. The estimated 
attach vector is the difference between the estimated output, 
y{t), and the actual, measured output, y{t). It is seen that the 
estimated attack vectors except the second one can converge to 
real attack vectors under no unknown inputs. However, when 
there do exist unknown inputs, the estimation of attack vectors 
is poor. This is because the CKF is not designed to inherently 
deal with unknown input vectors with unknown distributions. 


Fig. 5. Comparison between the actual states and estimated ones for the 
SMIB system utilizing the designed observer. 


C. Results for Observers 


The observer presented in Section V-B is presented concep¬ 
tually in Section VII In this section, we briefly present the 
solution of the observer gain in Section |V-C The Lipschitz 
parameters is set to unity, i.e., p = (p = y, = 1. Using CVX, 
a package for solving convex optimization problems ll33l on 
MATLAB, variables ei,e 2 and cr and a symmetric positive 
semi-definite matrix P, are computed; 


Cl = 0.0122, £2 = 0.0144, cr = 6.424, 
■ 0.4894 -0.0178 0.0627 -0.4632' 
-0.0178 0.0053 -0.0007 0.0067 
0.0627 -0.0007 0.7708 0.0274 

-0.4632 0.0067 0.0274 0.4945 


Then, the observer gain-matrix L is computed: 


L=fp- 




■ -6.023 15.9308 31.8673 12.0481 
-15.7483 42.5073 85.0299 31.5023 
4.2075 0.063 0.1261 -8.4166 

. -3.1177 8.6955 17.3942 6.2365 
a feasible solution for the linear matrix inequality O- 


The estimation of the states using the designed observer 
is shown in Fig. The observer succeeds in tracking the 
actual state under the presence of output attack vector and 


unknown inputs. The designed observer does not necessitate 


the knowledge of the initial conditions of the plant. In fact, 
other numerical simulations show that any observer initial 
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Fig. 6. Attack-vector estimation for the SMIB system utilizing the observer 
under unknown inputs and attack-vectors. 


conditions would still lead to accurate state estimation. 

Fig- HI shows the estimation of the predehned attack vector, 
v{t). As shown in Fig. the attack-vector estimator spectac¬ 
ularly accomplishes near-perfect disturbance estimation. Note 
that even though V 2 {t) is discontinuous, the estimator perfectly 
tracks the attack-vector value. 

Remark 7. The attack-vector estimator can be utilized to iden¬ 
tify the type of attack vectors and measurements’ disturbances 
against PMU outputs. Once estimated, the effect of cyber¬ 
attacks can be nullified by designing inputs that negate the 
effect of attack-vectors on state estimation. Furthermore, PMU 
channels experiencing high output disturbances or attacks 
might be diagnosed for potential malfunctioning or attacks, 
as illustrated by NESCOR failure scenarios (see Remark^. 

VII. Remarks and Insights 

Various functionalities of DSE methods and their strengths 
and weaknesses relative to each functionality are presented in 
this section based on (a) the technical, theoretical capabilities 
and (b) experimental results of the aforementioned techniques. 
Table represents a summary of this section. 

• Nonlinearities in the Dynamics: The developed estimation 
methods in this paper all operate on the nonlinear dynamics 
unlike the EKF that, by design, assumes a linearized version 
of the system dynamics. The presented observers in Sec¬ 
tion lY] assume, nonetheless, that the nonlinear component 
of the system dynamics is bounded, which is not a restricting 
assumption for most power systems. 

• Solution Feasibility: The main principle governing the de¬ 
sign of most observers is based on hnding a matrix gain(s) 
satisfying a certain condition — a solution to a matrix 
inequality or to a feasibility problem. Moreover, the gain- 
matrices for the two observers in Section |V] are solved via 
formulating a feasibility problem. The state estimates are 
guaranteed to converge to the actual ones if a solution to 
the LMI exists. In contrast, KE methods do not require that. 

• Unknown Initial Conditions: Unlike Kalman biters, most 
observer designs are independent on the knowledge of the 
initial conditions of the system; if the estimator’s initial 
condition is chosen to be reasonably different from the 


unknown, actual one, estimates generated from KF methods 
might not converge to the actual ones. For example, if the 
initial estimated state is chosen as a;^(0) = [3 Wo 0 0] , 
the CKF, UKF, and EKF estimates diverge. 

• Dynamics Subject to Unknown Inputs: As illustrated in 
the previous section, CKF-generated estimates is inferior in 
comparison with the estimation from the presented observer. 
Many derivatives of the KF are not designed to deal with 
varying unknown inputs, as the ones assumed in the paper. 
Similar conclusion applies for other KF methods. 

• PMU Measurements Subject to Cyber-Attacks: The observer 
for the nonlinear system outperforms the CKF, UKF, and 
EKF in the state estimation under attack-vectors as well as 
the estimation of the attack-vector. 


• Tolerance to Process and Measurement Noise: The ob¬ 
servers under consideration are tolerant to measurement and 
process noise similar to the ones assumed for KFs. By 
design, KF techniques are developed to deal with such noise. 

• Estimation Convergence Guarantees: Observers provide 
theoretical guarantees for the estimation convergence. How¬ 
ever, for KF there is no strict proof to guarantee that the 
estimation converges to actual states. 

• Numerical Stability: Observers do not have numerical stabil¬ 
ity problems while UKF can encounter numerical instability 
because the estimation error covariance matrix is not always 
guaranteed to be positive semi-debnite. 

• Tolerance to Parametric Uncertainty: KF-based methods 
can tolerate inaccurate parameters to some extent. Dynamic 
observers deal with parametric uncertainty in the sense that 
all uncertainties can be augmented to the unknown input 
component in the state dynamics (Ddif)). After redebning 
Dd{t) due to uncertainty, the observer is redesigned and 
the observer gain is updated to rebect the changes. 


Computational Complexity: The observers presented in the 
paper, CKF, UKF, and EKF all have computational com¬ 
plexity of 0{n^) IITaII . i20l . However, observers require 
fewer computational steps. Since the observers’ matrix gains 
are obtained offline by solving LMIs, observers are easier 
to implement as only the dynamics are needed in the 
simulations (e.g.. Equations (|^-(0 and ([T 2 I 1 ). Table [n| shows 
a comparison of the running time of CKE, UKE, and the 
observer from Section V-C| for the SMIB system — the 
estimation with the observer being the fastest among the 
three methods. 


TABLE II 

Running Time of the Three Estimators for SMIB System 



Running 

Time (seconds) 



CKF 

UKF 

Obs.2 (Sec. 

v-c[ 

a;(0) = 

[0 LUO 0 0]T 

3.28 

6.30 

0.76 


*'(0) = 

[3 LUO 0 0]^ 

3.27 

6.27 

0.76 



VIII. Conclusion 




























TABLE I 

Comparison between DSE techniques’ functionalities and characteristics (‘—’ indicates a lack of yes/no answer) 


F unctionality/Characteristic 

Kalman Filter Derivatives 

Robust Observers 

EKE 

UKF 

CKF (Sec. IV-C| 

Obs. 1 (Sec. V-B| 

Obs. 2 (Sec.|V-C| 

System’s Nonlinearities 

X 

/ 

/ 

✓ 

/ 

Feasibility 

— 

— 

— 

X 

X 

Tolerance to Different Initial Conditions 

X 

X 

X 

✓ 

✓ 

Tolerance to Unknown Inputs 

X 

X 

X 

✓ 

✓ 

Tolerance to Cyber-Attacks 

X 

X 

X 

X 

/ 

Tolerance to Process and Measurement Noise 

✓ 

/ 

/ 

✓ 

/ 

Convergence 

— 

— 

— 

✓ 

/ 

Numerical Stability 

— 

— 

— 

✓ 

/ 

Computational Complexity 

0(^3) 

0(n^) 

0(^3) 

0(n^) 

C>(n3) 


In this paper, we discuss different dynamic state estimation 
methods by presenting an overview of state-of-the-art estima¬ 
tion techniques and developing alternatives, including the cu- 
bature Kalman filter and dynamic observers, to address major 
limitations of existing methods. The proposed methods are 
examined on a single-machine infinite-bus benchmark system, 
under unknown inputs and attack-vectors. Based on both the 
technical, theoretical capabilities and the experimental results, 
we summarize the strengths and weaknesses of different 
techniques, and also provide the reader with recommendations 
for optimal DSE in power networks. 
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